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Abstract 

We show that two-photon transport is strongly correlated in one-dimensional waveguide cou- 
pled to a two-level system. The exact S-matrix is constructed using a generalized Bethe-Ansatz 
technique. We show that the scattering eigenstates of this system include a two-photon bound 
state that passes through the two-level system as a composite single particle. Also, the two-level 
system can induce effective attractive or repulsive interactions in space for photons. This general 
procedure can be applied to the Anderson model as well. 
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Creating strong photon-photon interaction at few-photon level is of great interest for 
quantum information sciences. In atomic gases, such interaction can be accomplished either 
with systems exhibiting electromagnetically induced transparency (EIT)^, or by reaching 
the strong-coupling regime of a two-level atom in a high-Q cavity^. However, in an on- 
chip, solid-state environment, which is crucial for practical applications, there have been 
significant challenges in implementing these concepts. For example, it is difficult to create 
the long-lifetime dark state, which is required for EIT effects, in most practical solid-state 
environments^. While the strong-coupling regime has been reached by placing a quantum dot 
in a high-Q photonic crystal microcavity^, doing so requires very accurate tuning of both 
the electronic and optical resonances to ensure simultaneous spectral and spatial overlaps^. 

In this Letter we propose and analyze in details an alternative scheme to create strong 
photon-photon interaction. Our approach exploits a unique one-dimensional feature for 
photon states in many nano-photonic structures. In a photonic crystal with a complete 
photonic band gap, for example, a line-defect waveguide forms a true one-dimensional con- 
tinuum for photons, since there is no other states within the gap. Here we show that by 
coupling a two-level system to such continuum, strong photon-photon interactions can be 
created (Fig. [1]). (Below we refer to the two- level system as the "atom"). In this system, the 
strong interaction arises from the fact that in a one- dimensional system, the re-emitted and 
scattered waves from the atom inevitably interfere with the incident waves. Moreover, since 
the atom, intuitively speaking, can at most absorb only one-photon at a time, the transport 
properties of multi-photon are strongly correlated. 

Compared with previous solid-state approaches, our scheme does not require the presence 
of long-lifetime dark state. Neither does this scheme necessitate detailed spectral tuning or 
spatial control of the two-level system, since it operates in the weak-coupling regime, and 
thus the one-dimensional continuum can be broadband. Moreover, the Hamiltonian of the 
system actually describes an exact photonic analogue of the Kondo effect, which is important 
for processing electronic quantum bits^. Our approach may therefore open a new avenue 
towards practical photon-based quantum information processing on-chip. 
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The system in Fig. [T] is modeled by the Hamiltonian^^: 

/( d d 

dx I -iVgC ] R (x)—c R (x) + iv g c\(x)—c L (x) 

+ VS( X ) (4(x)«7- + c R (x)a, + 4 (,),_ + 
+ E e a\a e + E g a g a g (1) 

where w g is the group velocity of the photons, and c\ i (x)(c\ j (x)) is a bosonic operator creating 
a right-going(left-going) photon at x. V is the coupling constant, a g (al) is the creation 
operator of the ground (excited) state of the atom, a + = a\a g (a_ = a g a e ) is the atomic 
raising (lowering) ladder operator satisfying cr+\n, — ) = \n, +) and a + \n,+) = 0, where 
\n, ±) describes the state of the system with n photons and the atom in the excited (+) or 
ground (— ) state. E e — E g (= Q) is the transition energy. This Hamiltonian describes the 
situation where the propagating photons can run in both directions, and is referred to as 
"two-mode" model. 

By employing the following transformation c\(x) = -4j(c|j(x) + c^(— x)), c\(x) = 
■^(c^(x) — c\{— x)), the original Hamiltonian is transformed into two decoupled "one-mode" 
Hamiltonians, i.e., H = H e + H , where 

H e = J dx(—i)v g cl(x) —c e (x) + J dxV5(x) (c\(x)<j- + c e (x)a + ) + E e a\a e + E g a g a g 

H Q = dx(-i)v g cl(x)—c (x). (2) 

H Q is an interaction-free one-mode Hamiltonian, while H e describes a non-trivial one-mode 
interacting model with coupling strength V = \[2V . H e is identical in form to the s-d 
modelii^, which describes the S-wave scattering of electrons off a magnetic impurity in 
three dimensions. Here, however, instead of a fermionic operators describing electrons, we 
have bosonic operators describing photons. 

The one-photon eigenstate for H e takes the form \k) = j dx[e lhx (8(—x) + tk6(x)) c*(x) + 
efccr + ]|0, — where 

is the transmission amplitude of magnitude 1, and = k _^ ir / 2 ^ s ^ ne excitation amplitude. 
The single photon experiences resonance when its energy k is close to the transition energy 
fl of the atom. For notational simplicity, v g and h are set to 1, and the subscript "e" in c\ 
is dropped hereafter. 



In this Letter we focus on the transport properties of the interacting Hamiltonian H e with 
two incident photons. For this Hamiltonian, as well as the Anderson model and the inter- 
acting resonance level model in condensed matter physics, attempts to diagonalize using the 
Bethe Ansatz have been published^^kL^. As we will emphasize below, however, a complete 
and correct description of the transport properties requires a careful re-examination of these 
solutions. In particular, the Bethe Ansatz solution constructed following the procedures in 

n 

Ref. [14j is in fact not complete for this purpose. Rather, to construct the scattering matrix, 
one needs one additional two-photon bound state. These can all be derived by a systematic 
approach detailed below. 

We first describe the general features of the scattering problem. Before and after the 
scattering, the photons are away from the atom, and thus the two-photon Hilbert spaces of 
the "m" (before scattering) and "out" (after scattering) states^ are the same space of free 
photons and consists of all symmetric functions of the coordinates of the photons, X\ and 
X2- This Hilbert space is spanned by a complete basis {\Sk, p ) : k < p} defined as 

(xi,x 2 \S kiP ) = ( e ^e^ 2 + e ikx *e ipx i) = ^e iEx ° cos (Ax) , (4) 

where E = k + p is the total energy of the photon pair, x c = \j2[x\ + X2), x = X\ — x 2 , 
and A = (k — E/2) = l/2(k — p) < 0. Alternatively, the same Hilbert space can instead be 
spanned by another basis : k < p] defined as 

(x 1 ,x 2 \A kp ) = — -^sgn(x) (e** 1 ^ 3 - e ikX2 e ipxi ) = — sgn(x) e lExc sin ( Ax) (5) 
2% y/2 2ix 

where sgn(x) = 8(x)—8(—x) is the sign function. We emphasize that, while both {|5jfc P ) : k < 
p} and {|Ajfc jP ) : k < p} are complete^, arbitrary linear combination {ak, p \Sk, p ) + bk )P \Ak tP ) : 
k < p] may not be. 

The transport properties of two photons, in the presence of the atom, are described by 
the S-matrix (S) that maps between the Hilbert space of the in and out states: |out) = S|in). 
The matrix element of the S-matrix, for example, (Sfc^lSISVy) is the transition amplitude 
of the process^. 

The S-matrix of the two-photon CclSGj clS will be derived below, can be diagonalized as 

S = ^t fc t P |W fe ,p}(W feiP | + ^t E \B E )(B E \, (6) 

k<p E 
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with 

^ ~ v / (fc _p )2 + r2 K* - P)l 5 **> + lT \ A ^ > 

Now we prove Eqs. (JSD-C3) by first showing that |Wfc P ) and \Be) are eigenstates of the 
scattering matrix. A two-photon eigenstate for H e has the general form: 

|$) = I / dx\dx 2 g(x\, x 2 )c (xi)<y (x 2 ) + / dx e{x)c\x)a + I |0, — ), (8) 



where e(x) is the probability amplitude of the atom in the excited state. Due to the boson 
statistics, the wavefunction satisfies g(xi,x 2 ) = +g(x 2 ,Xi). (g(xi,x 2 ) is continuous on the 
line Xi = x 2 for bosons.) 

From H e \<&) = E\<f>), we obtain the equations of motion: 

d d \ V 
~ % ~dx % ~dx E ) 9 ^ Xl ' X2 ' + ~2 ( e ( Xl ^( X2 ) + e { x 2)H x V) = °> 

-i^-E + ^ e(x) + V (g(0, x) + g(x, 0)) = 0, (9) 

where g(0,x) = g(x,0) e 1/2 x (g(Q~,x) + g(0 + ,x)). The functions g(xi,x 2 ) and e(x) are 
piecewise continuous. 

The interactions occur on the coordinate axes: x\ = 0, and x 2 = 0. Applying the 
equations of motions gives the following boundary conditions on the boundary of quadrants 
II (xi < < x 2 ) and III (xi, x 2 < 0): 

V 

-i [gfa, + ) - g(x u )) + —e{x x ) = 0, 

f -i^- - (£ - fi)) e(x0 + V(g( Xl , 0+) + ^(xi, 0")) = 0, (10) 
and on the boundary of quadrants II (x\ < < x 2 ) and I (0 < x±, x 2 ): 
-i (0(0+ x 2 ) - ^(0~, x 2 )) + ye(x 2 ) = 0, 

( -i^- - (£ - fi)) e(x 2 ) + \%(0+ x 2 ) + 0(0-32)) =0. (11) 
These boundary conditions must be supplemented by a further condition 

e(0-) = e(0+), (12) 
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which ensures the self-consistency. 

By boson symmetry we only need to consider the half space x\ < x 2 . In this half space, 
suppose g(xi,x 2 ) = B 3 e lkxi+ipX2 + A 3 e ipxi+,lkX2 for x\ < x 2 < 0, using Eqs. (jl0til2j) . we 
obtain g(x 1} x 2 ) = t k t p (B 3 e ikxi+ipX2 + A 3 e ipXl+ikX2 ) for < x x < x 2 , provided B 3 /A 3 = 
(k — p — iT)/{k — p + iT) as required from the continuity condition of e(x). Therefore, in 
the full quadrant III, the in-state, |Wjfc )P ) as defined by 

(x 1 ,x 2 \W k)P ) = (A 3 e ikxi+ipX2 + B 3 e ipxi+ikX2 ) B{x 1 - x 2 ) + (B 3 e ikxi+ipX2 + A 3 e ipxi+ikX2 ) 6{x 2 - x x ) 
oc (k - p)(x 1 ,x 2 \S k , p ) + ir(xi,x 2 \Ak,p), (13) 

is an eigenstate of the S-matrix with eigenvalue t k t p . This construction and the form of the 
solution is in essence the Bethe Ansatz methocU^ 1 ^. 

The set {|Wfe !P ) : k < p} however does not form a complete set of basis of the free 
two-photon Hilbert space. Instead, there exists one additional eigenstate of S-matrix, \Be), 
defined by Eq. (J7j). To see that \B E ) is an eigenstate of the S-matrix, suppose g{x\,x 2 ) = 
e iExc e -r\x\/2 - n q Uac } rari t ni, again using Eqs. fjlOH12p . we obtain g(xx,x 2 ) = tEe tExc e~ r ^ x ^ 2 
in quadrant I. Such bound state is important when calculating the ground-state energy in 
the Anderson modet^. We show here that it is also crucial to the scattering and transport 
properties. 

The set of eigenstates {jH^p, \Be)} forms a complete and orthonormal basis that 
spans the free two-photon Hilbert space. The orthonormality check is straightforward: 
(W k ,, p ,\W Kp ) = 5(k - k')5(p - p') = 6(A - A')5(E - E'), (B E \B E ) = 6(E - E'), and 
(Wk, p \BE) = 0. The completeness can be proven by checking that 

W = \W k , p )(W k , p \ + \Be)(B e \, (14) 

k<p E 

is indeed an identity operator. This, together with the eigenvalues t k t p and t E , prove Eq. 

We note that the two-photon bound state described by \Be), of which the spatial extent 
is 1/r, behaves as an effective single composite particle with an energy k + p, and remains 
integral when passing through the atom. The two-level system therefore provides the ca- 
pability of manipulating composite particles of photons^ without destroying them. This 
capability is important in quantum cryptography^ and quantum lithography^. 

For an arbitrary in-state of |in) = \S kliPl ), the momenta distribution of the out-state 
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(S k2 , P2 \out) is 



(S k2 , P2 \S\S kuPl ) = t kl t Pl S(A 1 -A 2 )5(E 1 -E 2 )+t kl t Pl 5(A 1 +A 2 )6(E 1 -E 2 )+B5(E 1 -E 2 ) (15) 



where the first two terms are the direct and exchange terms of each individual incident 
momentum; the third term with 



represents the background fluorescence due to the scattering. When Ai ^ A 2 , 
\B(E\, Ai, A 2 )| 2 is the probability density for the outgoing photon pair in (E 1 ,A 2 ) state, 
when the incoming photon pair is in (E\, Ai) state. 

The emergence of the background fluorescence is completely different from the well-known 
resonance fluorescence phenomenon where a strong laser beam is scattering off an ensemble 
of two-level systems^. In the current two-photon case, the background fluorescence results 
from the fact that the momentum of each photon is not conserved. Consequently the interac- 
tions with the two-level system redistribute the momenta of the photons over a continuous 
range, under the total energy and momentum conservation constraint. Furthermore, the 
locations of the poles in B, at k\ )2 = p\ )2 = Q — iT/2, correspond approximately to either 
one of the photons having an energy at Q. Thus, the background fluorescence arises as one 
photon inelastically scatters off a composite transient object formed by the atom absorbing 
the other photon. 

Fig. [2]plots normalized \B(E, Ai, A 2 )\ 2 for various photon pair energy E. \B(E, Ai, A 2 )\ 2 
is an even function of E — 2Q. When \E — 2Q\ < Y, there is a single peak centered 
at Ai = A 2 = 0. The height of the peak reaches maximum at E — 20, (Fig. 0(a)), 
and gradually decreases as \E — 2Q\ increases. When \E — 2Q\ = Y, the top of the 
peak becomes flat (Fig. (21(b)). When \E — 2Q\ > Y, there are four peaks centered at 
(±y/(E - 2VL) 2 - r 2 /2, ±^{E - 2VL) 2 - r 2 /2), respectively (Fig. (21(c) and (d)). For any E 
and Ai, the locations of the peaks for \B(E, Ai, A 2 )\ 2 are independent of Ai. In contrast, 
the 5-functions in the S-matrix (Eq. (fT5l) ) are located on the Ax = A 2 line. 

The emergence of the background fluorescence also manifests as an effective spatial in- 
teraction between the photons. For an in-state |in) = IS^Ai), the out-state is 



B 



lQiT 2 E 1 -2Vt + iV 



(16) 



7T [4A 2 - (Ei -2VL + 1Y) 2 ] [4A 2 - (E 1 - 2Vt + 1Y) 2 } ' 



(x c , x|out) = e 



,iE\x, 



,V2 



( 



t kl t Pl cos (Aix) 



4r 2 



,i(E 1 -2Q.)\x\/2-T\x\/2 



) 



4A 2 - 



(Ex -2VL + iY) 



2 



(17) 
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which takes the form e lElXc (x\(j)) , where (x\4>) is the wavefunction in the relative coordinate 
x. The deviation of the out-state wavefunctions from that of interaction-free case is large 
when Ai ~ ±(Ei/2 — Q), i.e., when at least one of the incident photons is close to resonance. 
Fig. E](a) plots the normalized deviation of \(x = O|0)| 2 from the interaction- free CclSG clS du 
function of E\ and A x . A positive (negative) deviation implies that the two photons bunch 
(anti-bunch) after scattering. The hyperbola 4 A 2 — {E\ — 2Q) 2 = T 2 indicate where the 
deviation is zero, thereby separate the bunching and anti-bunching regions. The deviation 
reaches maximum at E\ — 2Q = Ai = 0, when both incident photons are on resonance 
with the atom. The wavefunction for this case is shown in Fig. EJ^b), which exhibits the 
exponentially decaying feature in x. The two photons form a bound state after scattering, 
with half-width in space about 1/T. When E\ — 2Q is kept at zero, the height of the 
peak at x = decreases with increasing |Ai| (Fig. Efc)(d)). Fig. [3]^d) shows the case for 
Ai = —\/3T/2, where the peak at x = is completely depleted. Both bunching and 
anti-bunching behavior occur at other non-resonant E\ and Ai as well but is generally 
weaker. The resonance thus can induce either an effective repulsion or attraction between 
two photons. 

This work is supported by a Packard Fellowship. 
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FIG. 1: Schematics of the system. A two-level system is coupled to a one-dimensional continuum 
in which the photons, shown as wiggly waves, propagate in each direction. 
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FIG. 2: Background fluorescence as a function of A\ and A2 at various energy, (a) E = 0. (b) 
E = 2. (c) E = 4. (d) E = 6. i? = (r/2)J3, E = (E — 2fi)/(T/2), and A = A/(T/2). For any 
given the in- and out-states can be completely specified by one quadrant in the A1-A2 plane. 
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FIG. 3: (a) Plot of the normalized deviation of the relative wavefunction: | (x = 0|</>)| 2 / (\/2/27r) 2 — 1. 
x = (T/2)x. The two gray lines indicates the zero value. (b)(c)(d): plots of |(x|^>)| 2 for E\ = 0. 
(b) Ai = 0. (c) Ai = —0.8. (d) Ai = —y/3. Gray lines indicate the interaction-free case. 
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